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Abstract 

A theoretical and experimental Investigation 
of single-phase and particle- laden weakly swirling 
jets was conducted. The jets were Injected ver- 
tically downward from a 19 mm diameter tube with 
swirl numbers ranging from 0 to 0.33. The 
particle-laden jets had a single loading ratio 
(0.2) with particles having an SMD of 39 ym. Mean 
and fluctuating properties of both phases were 
measured using nonlntruslve laser based methods 
while particle mass flux was measured using an 
Isokinetic sampling probe. The continuous phase 
was analyzed using both a baseline k-c turbulence 
model and an extended version with modifications 
based on the flux Richardson number to account for 
effects of streamline curvature. To highlight 
effects of Interphase transport rates and particle/ 
turbulence interactions, effects of the particles 
were analyzed In three ways, as follows: 

(1) locally homogeneous flow (LHF) analysis, 
where Interphase transport rates are assumed to 
be Infinitely fast; (2) deterministic separated 
flow (DSF) analysis, where finite Interphase 
transport rates are considered but particle/ 
turbulence Interactions are Ignored; and (3) sto- 
chastic separated flow (SSF) analysis, where both 
effects are considered using random walk 
computations . 

Nomenc 1 at u re 

a acceleration of gravity 

C 0 drag coefficient 

parameters In turbulence model 

d Injector diameter 

d particle diameter 

P 

f particle mass fraction 

G particle mass flux 

k turbulence kinetic energy 
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L e dissipation length scale 

nip particle mass 

n number of particle groups 

n^ number of particles per unit time In 

group 1 

p pressure 

r radial distance 

R Injection tube radius 

Re Reynolds number 

S source term 

S source term due to particles 

P4> 

t time 

t eddy lifetime 

u axial velocity 

v radial velocity 

Vj ’volume of computational cell j 
w angular velocity 

x axial distance 

x particle position vector 

AXp relative path length of particles In an eddy 

At time of particle residence In an eddy 

c rate of dissipation of turbulence kinetic 
energy 

y^ turbulent viscosity 

P density 

turbulent Prandtl/Schmldt number 
4 > generic property 
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Subscripts 

c centerline quantity 

p particle property 

o Injector exit condition 

® ambient condition 

m maximum value 

Supersc r ipts 

{ )' root-mean square fluctuating quantity 

(“) time averaged value 

(*) vector quantity 

(~) Favre averaged value 

Intr oduction 

Multiphase flows are extremely Important In 
many practical devices. In a gas turbine engine, 
for example, fuel Is sprayed Into a highly tur- 
bulent, recirculating flowfleld, where It evapo- 
rates and takes part In chemical reactions. 
Particle-laden flows are being studied as a step 
toward a better understanding of complex multi- 
phase flows such as those In future advanced air- 
breathing engines. Particle-laden flows allow the 
study of Interactions between the continuous and 
dispersed phases without Interference from vapor 
Izatlon or combustion effects. The objective of 
this Investigation was to extend previous work on 
free particle-laden jets to consider effects of 
swirl, which Is an Important aspect of many prac- 
tical sprays. The performance of typical two- 
phase flow models of this problem was of 
particular Interest. 

Single-phase swirling jets have been studied 
several times In the past. 1 "* These flows are 
normally characterized by their swirl number. For 
a free jet In stagnant surroundings, the swirl 
number can be expressed as follows: 

2 

/ puwr dr 

s « — — - 1 * 2 ( i ) 

R f 0 jpu * (p - pjj r dr 

where S Is an Invariant of these flows. 5 For 
swirl numbers greater than roughly 0.6, the 
adverse pressure gradient set up by the decay of 
angular velocity causes reversal of the streamwlse 
velocity and formation of a zone of recirculation 
along the axis. However, even at lower swirl 
numbers, the presence of swirl has large effects 
on the structure of a jet. For example. 

Increasing swirl numbers cause Increased rates of 
jet growth, entrainment of ambient fluid, and 
decay of streamwlse velocity. 

Earlier theoretical and experimental Inves- 
tigations of particle-laden jets have recently 
been reviewed; 6 therefore the present discussion 
of past work will be brief. The present study Is 
an extension of earlier Investigations 
Involving analysis and experiments for particle- 
laden jets. Three methods of analyzing multiphase 


flows were considered, as follows: (1) a locally 

homogeneous flow ( LHF) model, where properties of 
both phases were taken to be the same; (2) a 
deterministic separated flow (OSF) model, where 
finite Interphase transport was considered but 
effects of turbulence on particle motion were 
Ignored; and (3) a stochastic separated flow (SSF) 
model where effects of both finite Interphase 
transport rates and turbulence on particle motion 
were considered using random- samp! Ing techniques. 
In general, the LHF and DSF models over- and 
underestimated rates of particle spread and flow 
development, respectively. In contrast, the SSF 
model yielded encouraging predictions of flow 
structure- except at high particle mass loadings, 
where effects of particles on turbulence prop- 
erties (termed turbulence modulation by A1 Taweel 
and Landau 11 ) which were not considered In the 
theory, were felt to be responsible for the 
deficiency. 

In the present Investigation, the results of 
the previous studies were extended to particle- 
laden, weakly swirling jets. Relatively few 
theoretical studl'es treating particle- laden 
swirling flows have been published. Hamed 12 
modeled the trajectories of solid particles In a 
flow downstream of swirl vanes. Drlng and Suo 13 
calculated particle trajectories In a free-vortex 
flow. Effects of turbulence on the particle tra- 
jectories, however, were not considered In either 
of these studies. Furthermore, no previous 
experimental studies of particle- laden Jets with 
swirl could be found In the literature. 

Due to the absence of existing measurements, 
the present experimental study was conducted to 
obtain data that could be used to assess the pre- 
viously described models In a swirling flowfleld. 
In order to obtain significant Interaction between 
phases, relatively small solid particles (with an 
SMD of 39 vim) were used. Emphasis was placed on 
obtaining gas-phase velocities In the presence of 
particles, without signal Interference from the 
particles. Initial conditions were measured near 
the jet exit. Single-phase jets were also studied 
to establish baseline results, prior to treating 
the particle-laden flows. 

E x£e rj ment al M ethods 

Test Ap paratus 

Both the single-phase and particle-laden jets 
were directed vertically downward within a large 
(1.8 m-square x 2.4 m high) screened enclosure. 

The injector could be traversed vertically within 
the enclosure while the enclosure and Injector 
could be traversed together In the two horizontal 
directions. This arrangement allowed rigid 
mounting of all optical Instrumentation used 
during the study. 

The jet tube had an Inside diameter of 19 mm 
and extended vertically downward for 100 Injector 
diameters. Swirl was generated by Injecting air 
tangentially through four 9.5 mm long slots, 
located 74.6 mm upstream of the jet exit. The 
swirl number was changed by varying the amount of 
air Injected through the tangential slots. The 
air flowing through the injector tube was fil- 
tered. Both the swirl and main air streams were 
metered with calibrated, critical-flow orifices. 
The solid glass particles used during the study 
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were Injected Into the flow far upstream of the 
Injection tube using a vibrating, variable-speed 
screw feeder. Particle properties and test con- 
ditions are summarized In Table I. 

Instrumentation 

Gas velocity . Mean and fluctuating gas 
velocities In the single-phase jets were measured 
with a two-channel (frequency shifted on both 
channels) Ar-Ion laser veloclmetry (LV) system. 

The backscatter mode was used with a collecting/ 
receiving lens focal length of 750 mm. A beam 
expansion ratio of 3.75 was used to Improve 
slgnal-to-nol se ratio. Both the jet and the 
ambient surroundings Inside the enclosure were 
seeded (with 1 ym nominal diameter aluminum oxide 
particles) to eliminate concentration bias near 
the edge of the jet. Two horizontal traverses, 

90° apart, at each axial position yielded all 
required velocity components, since two velocity 
components were measured simultaneously. 

Gas-phase velocities were measured In the 
presence of particles using a single-component 
phase-doppler anemometer, described In detail by 
Bachalo and Houser. 14 With this Instrument, 
particle size and velocity were measured simul- 
taneously. Gas-phase velocity was measured by 
seeding the particle-laden jet and ambient sur- 
roundings with 1 ym nominal aluminum oxide par- 
ticles and using the velocity measured for this 
size range to represent the continuous phase. The 
green (514.5 nm) line of the Ar-Ion laser was used 
with a 602.4 mm focal-length lens. To reduce the 
size of the probe volume, a beam expansion ratio 
of 3 was used. The receiving optics were mounted 
30° off-axis In the forward scatter direction, 
with scattered light collected using a 495 mm 
focal-length lens. Since frequency shifting was 
not available for this system; measurement of 
radial and angular velocities were not performed, 
and only mean and fluctuating axial velocities are 
reported . 

Particle velocit y. Particle velocities were 
measured using the two- component LV system used 
for the single-phase flows, except that the 
receiving optics were placed off-axis approxi- 
mately 30° In the forward scatter direction. The 
collecting lens for the receiving optics had a 
focal length of 602.4 mm and frequency shifting 
was used for both channels. Present data are 
number averages over all particle sizes. During 
the particle velocity measurements, seeding par- 
ticles were not Introduced Into the flow and both 
laser power and detector gain were reduced to 
further Insure that only signals from the test 
particles were observed. At each axial position, 
two traverses, 90° apart, yielded all required 
velocity components. 

Particle mass flux . Particle mass flux was 
measured using an Isokinetic sampling probe. 
Samples were collected on filter paper for a timed 
Interval and weighed. Probes having Inside diam- 
eters of 2 and 5 mm were used to Insure adequate 
resolution and reasonable sampling times In var- 
ious regimes of the flow. Measured particle mass 
fluxes Integrated across the jet were within 
±10 percent of the calibrated particle flow rate 
at all axial locations. 


Theoretical Methods 
General Description 

The analysis Is limited to steady, axlsym- 
metrlc, dilute, sol Id-particle-laden, weakly 
swirling turbulent jets In an Infinite stagnant 
media. The swirl number, calculated from Eq. (1) 
Is restricted to values less than approximately 
0,5 to prevent any zones of recirculation. The 
boundary-layer approximations are adopted, how- 
ever, the radial pressure gradient, which Is 
usually neglected In the boundary- layer analysis. 
Is considered. A k ~c turbulence model Is used 
to provide closure since It has modest computa- 
tional requirements. Effects of streamline cur- 
vature on the k-e model are considered. The 
Injector exit Mach number Is less than 0.3; 
therefore the kinetic energy and viscous dissipa- 
tion of the mean flow are neglected with little 
error . 

Three methods of treating multiphase flow, 
typical of current practice, are considered. The 
methods are: (1) locally homogeneous flow (LHF), 

where Interphase transport rates are assumed to 
be Infinitely fast and the flow can be treated 
like a single-phase, variable-density fluid; 

(2) deterministic separated flow (DSF), where 
finite Interphase transport rates are considered 
but the dispersed phase Is assumed to Interact 
only with mean properties of the continuous phase; 
and (3) stochastic separated flow (SSF), where 
Interphase transport rates and effects of turbu- 
lent dispersion of the dispersed phase are 
treated. All three methods will be only briefly 
discussed since they have been fully described 
elsewhere. 

Continuous Phas e 

Mean quantities for the continuous phase are 
found by solution of governing equations for con- 
servation of mass and momentum In conjunction with 
second-order turbulence model equations for tur- 
bulent kinetic energy and Its rate of dissipation. 
The volume fraction of the particle phase was 
neglected, since void fractions for the present 
flows exceeded 99.8 percent. The governing equa- 
tions for the continuous phase can be put Into the 
following general form: 

fe <**> * f k < r ^> 



The parameters 4>, S^, and Sp^, appearing In 
Eq. (2), as well as empirical constants, are sum- 
marized In Table II. A modification to the c 
equation, based on use of a flux Richardson number 
correction to account for the effect of streamline 
curvature, 15 was also studied. For this modi- 
fication, C c ] Is replaced by 

C cl (1 + °* 9 R f> {3) 
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where 



A consequence of the presence of angular 
velocity Is that even though the standard boundary 
layer assumptions have been made, the radial 
momentum equation still must be considered: 

1 

— _ / c \ 

3r - r v ; 

and cross-stream pressure gradients are not neg- 
ligible. Because of the decay of angular veloc- 
ity, the axial pressure gradient Is also Included. 


The present approach Involves modification and 
extension of methods reported by Gosman and 
Ioannldes. 16 Properties are assumed to be uni- 
form within each eddy and to randomly change from 
one eddy to the next. At the start of a particle- 
eddy Interaction, the velocity of an eddy Is found 
by making a random selection from the probability 
density function (PDF) of velocity- assuming an 
Isotropic Gaussian PDF having standard deviations 
(2k/3)'/ 2 and mean values TJ, V, and w. A 
particle Is assumed to Interact with an eddy for 
a time which Is the minimum of either the eddy 
lifetime or the time required for a particle to 
cross an eddy. These times are estimated fol- 
lowing past practice, 7~10 assuming that the 
characteristic size of an eddy Is the dissipation 
length scale 


L =C 3 / V /2 /e 
e y 


(9) 


The turbulent viscosity was calculated as 
usual : 



The flow leaving the Injector was similar to 
fully developed flow and had no potential core. 
The boundary conditions for Eq. (2) are: 


and that the eddy lifetime Is 

1/2 

t e = L e /(2k/3) (10) 

Therefore, particles are assumed to Interact with 
an eddy as long as both the time and distance of 
Interaction satisfy the following criteria 


4 > = rw : r * 0, = 0; r -> «>, $ = 0 

a r 

(7) 

<$> = rw : r = 0, rw = 0 ; r -» », rw = 0 

Initial conditions were measured at 
x/d = 0.5. For the single- phase cases, all three 
mean and fluctuating velocities were measured. 
c 0 was calculated from the definition of a 
turbulent length scale, as follows: 

C k 3/2 

. - -X- /fli 


where L was chosen as a fraction of the Initial 
jet width to provide good agreement with the 
Initial axial profile for k. Since only U and 
u 1 could be measured for the part Icle- laden 
cases, Initial values of the continuous phase 
angular velocity were estimated by subtracting the 
measured particle angular momentum from the single- 
phase values. Single-phase Initial values of k 
and e were used for the particle- laden jet 
calculations . 

D|spe rsed Pha se_ ( SSJL Iprmul at Ion) 

The dispersed phase was treated by solving 
Lagrangian equations for the trajectories of a 
statistically significant sample of Individual 
particles (n groups defined by Initial position, 
velocity, and sample) as they move away from the 
injector and encounter a random distribution of 
turbulent eddies. This approach provides a means 
of treating effects of turbulent fluctuations on 
particle drag and dispersion as well as effects 
of particles on turbulence properties. 

Key elements of the SSF model are the methods 
used to specify eddy properties and the time of 
Interaction of a particle with a particular eddy. 


At < t , 

p - e' 


Ax„ 


< L 


(ID 


Assumptions for particle trajectory calcula- 
tions are typical of analysis of dilute particle- 
laden flows:® drag Is treated empirically, 
assuming quasisteady flow for spherical particles 
with no Influence of nearby particles; particle 
collisions are neglected; since p p /p > 200 for 
present tests, effects of virtual mass, Bassett 
forces and Magnus forces are neglected with little 
error; and static pressure gradients are negli- 
gible. Local ambient properties are fixed by 
Instantaneous eddy properties, as described 
earlier, which Implicitly provides for effects of 
turbulent fluctuations on particle dispersion and 
drag. 


With these assumptions, the position and 
velocity for each particle group can be found by 
Integrating 



where 1 = 1,2,3 and the velocities shown In these 
equations are Instantaneous velocities for a par- 
ticular eddy and particle group. 5 Since the 
particles used In the study were relatively round, 
the standard drag coefficient for solid spheres 
was taken as follows:** 

c 0 « re ( 3 ♦ -6“-)’ Re < 100 °; 

C D = 0.44 , Re > 1000 (14) 
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Particle Source Terms 


The Interaction between particles and the 
continuous phase yields source terms In the gov- 
erning equations for conservation of axial and 
angular momentum. The source terms are found by 
computing the net change In momentum as each par- 
ticle group 1 passes through computational 
cell j 


pu. 




(15) 



tl -l 


£ 


1=1 




(16) 


where n^ Is the number of particles per unit 
time In each group. 


an entrainment Increase of 2 percent. The angular 
and radial momentum equations required modifica- 
tions to the standard solution procedure. GENMIX 
utilizes an expanding grid In the streamwlse 
direction and downstream values of r are not 
known before the axial momentum equation Is 
solved. In order to preserve the marching char- 
acter of the solution procedure, the angular 
momentum equation Is solved first to obtain rw 
and then the radial momentum equation Is Inte- 
grated across the flowfleld using downstream 
values of rw and upstream values of r. Then 
ap/ax Is calculated and a correction applied to 
It based on conservation of axial momentum. This 
procedure has been previously reported by 
Slddhartha. 18 It was found to conserve axial 
momentum within 2 percent. The dispersed phase 
was computed using a second-order finite- 
difference algorithm for both separated flow 
models. For the SSF model, 2000 and 4000 particle 
groups were tracked through the flowfleld for the 
non-swirl and swirl cases, respectively. 


DSF and LHF Models 


Results and Discussion 


DSF model . Effects of turbulent fluctuations 
on particle drag and dispersion as well as effects 
of turbulence modulation are Ignored for the DSF 
model. Particle trajectories are found by Inte- 
grating Eqs . (10) and (11), but the local mean 
velocity of the continuous phase replaces the 
Instantaneous eddy velocity. Each Initial condi- 
tion yields a single deterministic trajectory; 
therefore, 1000 particle groups suffice to numer- 
ically close the solution. Effects of particle 
drag In the mean momentum equation are found from 
Eqs. (15) and (16), similar to the SSF 
calculations . 

LHF model . This approximation Implies that 
both phases have the same instantaneous velocity 
at each point In the flow; therefore, the flow 
corresponds to a variable-density single-phase 
fluid whose density changes due to changes In 
particle concentration. Turbulent dispersion of 
particles Is then equivalent to that of a gas and 
particle Inertia fully Influences turbulence 
properties; l.e., the method Implicitly accounts 
for effects of turbulence modulation to the extent 
that the no-slip assumption Is correct. 

The treatment of the variable-density fluid 
Is similar to past practice, & however f Is 
defined here as the mass fraction of particles In 
the fluid. Measured Initial values of f were 
used for the LHF predictions. Through the 
assumption of no-slip, there Is no need to compute 
particle trajectories and all particle source 
terms In the governing equations for the 
continuous phase are zero. 

Numerical Solution 

The calculations for the continuous phase 
were performed using a modified version of 
GENMIX. 18 The computational grid for the cases 
without swirl was similar to past work: 6 ' 10 
33 cross-stream grid nodes and streamwlse step 
size was limited to 5 percent of the current flow 
width or an entrainment Increase of 5 percent - 
whichever was smaller. For the cases with swirl, 
33 cross-stream grid nodes were also used but 
streamwlse step size was reduced to the smaller 
of either 2 percent of the current flow width or 


Single-Phase Results 

The single-phase results will only be briefly 
discussed since the present emphasis Is on 
particle-laden flows. Typical results are shown 
to give an Indication of the effect of swirl on 
the gas-phase flowfleld. Predictions are shown 
both Ignoring and considering the curvature cor- 
rection to the dissipation equation based on the 
flux Richardson number. 

Axial variation of flow pr opert ies . Predicted 
and measured mean streamwlse velocity along the 
axis, maximum angular velocities, and turbulent 
kinetic energy along the axis for the single-phase 
jets are Illustrated In Figs. 1(a) to (c). 

Results for swirl numbers of 0, 0.19, and 0.33 are 
shown on the figures. The rate of decay of 
streamwlse velocity along the axis Is Increased 
as swirl Is Increased. Turbulent kinetic energy 
also Increases quite dramatically with Increasing 
swirl. For the S = 0.33 case, turbulent kinetic 
energy rises steeply almost Immediately and peaks 
at an x/d of only 2 before decreasing once 
again. The maximum angular velocity decays quite 
rapidly, reaching one-half Its original value In 
less than 2 Injector diameters for S = 0.33. 
Predictions are In reasonably good agreement with 
the measurements at low swirl. Predictions are 
shown for both the standard k -c model and the 
flux Richardson correction. As shown In 
Fig. 1(b), the curvature correction greatly 
Increases k Immediately downstream of the 
Injector exit. However, for S = 0.33, the rapid 
rise of k near the Injector exit Is not pre- 
dicted very well, even with the curvature correc- 
tion. As shown In Figs. 1(a) and (c), the 
Increase In k using the curvature correction 
Increases the rate of decay of both axial and 
angular velocities. 

Radial profiles . Radial profiles of mean and 
turbulence quantities at x/d = 15 for the jet 
without swirl are Illustrated In Fig. 2. The jet 
width Is predicted reasonably well, however, cen- 
terline values of k and peak values of Reynolds 
stress are somewhat overestimated. Radial pro- 
files at x/d = 30 are not shown, however. 
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agreement between predictions and measurements was 
quite good. 

Radial profiles at x/d = 5 for the two 
single-phase jets with swirl are shown in 
Figs. 3(a) and (b). As expected. Increasing the 
swirl number Increases the rate of spread of the 
jet. Turbulent kinetic energy and Reynolds 
stress also Increase with Increasing swirl at this 
axial location. Predictions are reasonably good 
with the flux corrected version yielding best 
agreement with measurements. Radial profiles at 
x/d = 10 for the same swirling flows are Illus- 
trated In Figs. 4(a) and (b). Findings are quite 
similar to those obtained at x/d = 5 (Illustrated 
In Figs. 3(a) and (b)) . 

Partlcl e_- L aden Jets 

Predictions and measurements of the particle- 
laden Jets are discussed In the following. As 
discussed earlier, only mean and fluctuating axial 
velocities could be measured for the gas phase; 
therefore. It was necessary to estimate Initial 
values of angular velocity and turbulent kinetic 
energy for the gas phase. Initial angular veloc- 
ity of the continuous phase was estimated by sub- 
tracting the measured Initial particle- phase 
angular momentum from the values obtained for the 
single-phase cases. For both swirl flows, 
particle-phase angular momentum was approximately 
10 percent of the values measured for the single- 
phase cases. For the predictions shown here. 
Initial values of k were assumed to be the same 
as the single-phase flows. Measured values of 
u' for the gas phase were approximately 20 per- 
cent lower across the entire jet width for the 
particle- laden cases (at the Initial condition of 
x/d = 0.5) than for the corresponding single-phase 
flows. Predictions showed that reductions of k 
of 20 percent caused negligible changes In flow 
properties except very close to the Injector. 
Initial values of e for the particle laden jet 
predictions were also unchanged from the single- 
phase cases. Only the standard k-e model was 
utilized for the predictions reported In the 
following. 

Axi al vari ation of flow properties . Predicted 
and measured mean gas- phase axial velocities are 
shown In Fig. 5 for swirl numbers of 0, 0.16, and 
0.3. For no swirl, axial velocity decays more 
slowly for the parti c le- laden flows than for the 
single- phase flow due to the momentum exchange 
from the particles. For the swirl flows, inlet 
swirl numbers were reduced approximately 10 per- 
cent for the particle- laden cases; therefore, a 
direct comparison with the single-phase measure- 
ments cannot be made. Predictions from both the 
SSF and LHF models are nearly identical and show 
good agreement with the data for both swirl flows. 

Predictions and measurements of mean particle 
axial and maximum angular velocities as a function 
of streamwlse distance are Illustrated In Figs. 6 
and 7. As expected, the LHF model overestimates 
the rate of decay of velocity for all cases since 
the particles possess significant Inertia. Dif- 
ferences between the DSF and SSF models were small 
with both models giving reasonable agreement with 
the measurements. For both cases with swirl, 
particle velocity was lower than the gas phase at 
x/d = 0.5. Particle axial velocity Increases due 
to momentum exchange from the continuous phase 


before finally beginning to decay, c.f. Fig. 6. 

It Is encouraging that both separated- fl ow models 
correctly predict this behavior. 

Radial vari ation of flow pr operties . The 
radial variation" of both gasand particle-phase 
flow properties at x/d = 5 and 30 are Illustrated 
In Figs. 8 and 9. Similar to past results, 6 
gas-phase predictions are nearly Identical for 
both the SSF and LHF models since particle loading 
ratio Is relatively low. Predicted particle 
velocity shows the largest variation between 
models at x/d = 5. Here, the LHF model under- 
estimates particle velocity while both separated- 
flow models are In reasonably good agreement with 
the data. At x/d = 30, shown In Fig. 9, the 
different physical assumptions embodied In the 
models are quite evident In the plot of radial 
mass flux. The LHF model overestimates the dis- 
persion of particles while the DSF model clearly 
underestimates particle dispersion. Only the SSF 
model correctly predicts particle mass flux at 
this location. The predictions of fluctuating 
particle properties for the SSF model, also shown 
In Figs. 8 and 9 are reasonably good, however, 
u'p Is consistently underestimated. This Is 
probably caused by the assumption of Isotropic 
fluctuations In the SSF model and the monodlsperse 
assumption In the predictions, since the Injected 
particles have a standard deviation of 15 ym. 
Measured values of u* were always greater than 
v' and w 1 for the single-phase jet and It Is 
expected that this behavior should be similar for 
the particle- laden jet. 

Radial profiles of flow properties In the 
swirling, particle- laden Jets are Illustrated In 
Figs. 10 to 14. Gas phase properties at x/d = 5 
are shown In Fig. 10 for swirl numbers of 0.16 and 
0.3. Predictions using the SSF and LHF models are 
nearly identical and show that the width of the 
jet is slightly overestimated. Predictions this 
close to the Injector are somewhat dependent upon 
Initial conditions. Turbulence modulation by the 
particles, not Included In these predictions, 
could also be a factor In the decreased width of 
the jet. 

Radial profiles of particle properties at 
x/d = 5 are Illustrated In Figs. 11(a) and (b). 
Both separated- flow models slightly overestimate 
axial and angular particle velocities because the 
gas phase width Is overestimated. The LHF model 
underestimates particle velocities here since the 
particles have Inertia and do not Immediately 
respond to the continuous phase. Particle flux, 
also shown In Fig, 11, is quite different for the 
swirl cases than for the round jet without swirl. 
The swirl component has shifted the maximum flux 
outward from the center of the jet. Differences 
between the models are clearly evident In particle 
flux predictions. Even though initial values of 
particle concentration were used as Initial con- 
ditions for the LHF model, the high rate of tur- 
bulent dispersion due to the no-slip assumption 
has caused the predicted maximum particle flux to 
shift to the center of the jet. This Is clearly 
not correct. Both the SSF and DSF models predict 
an off-center maximum, however the predicted peak 
Is shifted radially outward from the measured 
value due to the overprediction of angular veloc- 
ity. The effect of turbulent dispersion Is 
clearly Important for this flow since the par- 
ticles have spread much farther than the DSF model 
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predicts. Fluctuating particle velocities at 
x/d = 5 are also Illustrated In Figs. 11(a) and 
(b) and show distinct differences for the two 
swirl numbers. For the lower swirl number, 

S = 0.16, SSF predictions for all three fluc- 
tuating velocities compare quite well with data. 
For the S = 0.3 case, however, u' p Is somewhat 
overestimated. For the single-phase case at 
x/d = 5, the measured peak value of u' Is 
approximately 40 percent higher than measured for 
the particle-laden jet. Some of the reduction In 
u 1 Is due to the reduction In swirl number caused 
by the Introduction of particles (approx. 10 per- 
cent), however part of the reduction Is thought 
to be due to turbulence modulation by the par- 
ticles. This appears to be more Important for the 
swirl flows. 

Radial profiles of gas-phase properties at 
x/d = 20 are Illustrated In Fig. 12 for both swirl 
flows. Predictions are In better agreement with 
measurements at this position than closer to the 
Injector, although the jet width Is still slightly 
overestimated. Predictions at x/d = 20 are not 
as sensitive to Initial conditions and the swirl 
component has almost completely decayed. Again, 
there Is little difference between predictions of 
gas-phase properties for the LHF and SSF models. 

Radial profiles of particle properties for 
both swirl flows are Illustrated In Figs. 13(a) 
and (b). Predictions of axial velocity are In 
good agreement with measurements for all three 
models. Predicted values of u r p underestimate 
the measurements while v‘ p and w‘ p are In 
reasonably good agreement with measurements using 
the SSF model. Since effects of swirl have 
decayed at this axial location. Ignoring the ani- 
sotropy of the continuous phase and using a single 
particle size for the predictions are the main 
reasons for this behavior- as discussed earlier. 
The particle mass flux predictions again highlight 
the different physical assumptions of the three 
models. Particle mass flux measurements Indicated 
that between x/d = 5 and x/d = 10, the maximum 
mass flux shifted to the center of the Jet for 
both swirl cases. Since angular and radial 
velocities have decayed to relatively small values 
at this distance and would tend to move particles 
outward, the only mechanism for transport Inward 
Is turbulent dispersion. As shown In Fig. 13(a), 
for the lower swirl number, the SSF model predicts 
this shift In mass flux. In contrast, the DSF 
model predicts a very narrow distribution with 
almost no particles at the center of the jet. For 
S = 0.3, the predicted maximum particle flux for 
the SSF model has not completely shifted to the 
center, however It Is clearly evolving In this 
direction. Again, the DSF model predictions show 
the particles remaining In a narrow area. In 
contrast to the no-swirl case, the LHF model 
underestimates particle dispersion for both swirl 
flows at x/d = 20. This behavior Is caused by 
neglecting the angular Inertia of the particles 
which tends to transport them radially. 


Conclusions 

Major conclusions concerning the models that 
were evaluated during this Investigation are as 
follows : 

1. For single-phase weakly swirling jets, a 
correction to the dissipation equation, based on 
the flux Richardson number and designed to account 
for the effect of streamline curvature, gave bet- 
ter agreement with measurements than the standard 
k-e turbulence model. This modification always 
Increased predicted levels of k across the 
entire width of the jet which Increased both the 
rate of spread of the jet and the rate of velocity 
decay with streamwlse distance. 

2. The SSF model, which accounts for both 
particle Inertia and particle/turbulent Interac- 
tions, yielded reasonably good results for the 
particle-laden jets. The prescription of eddy 
properties was not changed for these calculations 
from Its original calibration. 8 The DSF model 
performed quite poorly for prediction of the par- 
ticle mass flux distribution for the weakly 
swirling jets considered In this study. This 
suggests that turbulent dispersion of the par- 
ticles was Important for the present flows. For 
the swirling particle-laden jets, In contrast to 
the particle-laden jets without swirl, the LHF 
model underestimated the spread of the dispersed 
phase. Thus, the angular and radial Inertia of 
the particles are Important considerations In the 
mass flux distribution of the particles for the 
present flows. 

3. For the swirling particle-laden Jets, the 
continuous phase jet width was somewhat overesti- 
mated using the standard k-c turbulence model. 
This Is In contrast to the single-phase flows, 
where Jet widths were underestimated with the 
standard k-c model. Turbulence modulation, not 
considered In the predictions presented here, was 
probably a significant factor In the reduction of 
the jet width. 
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TABLE I. - SUMMARY OF TEST CONDITIONS 3 


Parameter 

Single-phase jets 

Particle-laden jets 


Case 1 

Case 2 

Case 3 

Case 1 

Case 2 

Case 3 

Centerline air axial velocity, m/s b 

14.86 

12.94 

13.34 

13.75 

11.9 

12.2 

Maximum air angular velocity, m/s b 

0 

3.16 

5.92 

0 

— 

— 

Swirl number 0 

0 

.19 

.33 

0 

.16 


Centerline particle axial velocity, m/s b 

— 



13.8 

10.39 

10.26 

Maximum particle angular velocity, m/s b 




0 

1 .48 

2.26 

SMD, pm d 




39 

39 

39 

Loading rat1o e 




.2 

.2 

.2 


a Amb1ent temperature and pressures, 296. K, 97 kPa; Injector Inside diameter, 19 mm; 

particle density, 2500 kg/m 3 , 

^Measured at x/d = 0.5. 

Calculated from Eq. (1). 

d The size distribution has a standard deviation of 15 pm. 
e Rat1o of Injected particle mass flow rate to air mass flow rate. 
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(a) Streamwise velocities. 

Figure 1. - Variation of flow properties in the streamwise 
direction for single-phase jets. 




(b) Turbulent kinetic energy. 
Figure 1. - Continued. 







(c) Maximum angular velocities. 
Figure 1. - Concluded. 





Figure 2. - Radial profiles of flow properties in 
the single-phase nonswirling jet at x/d * 15. 
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Figure 5. - Streamwise gas velocities along the axis of the 
particle-laden jets. 



Figure 6. - Particle streamwise velocities along the axis of the 
particle-laden jets. 
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Figure 9. - Radial variation of flow properties in the nonswirling 
particle-laden jet at x/d = 30. 
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Figure 11, - Radial profiles of particle properties in 
swirling particle-laden jets at x/d = 5. 
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Figure 13. - Radial profiles of particle 
properties In swirling particle-laden 
jets at x/d = 20. 
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